Mem. S.A.It. Vol. 73, 23 

© SAIt 2002 Memorie della 

ROSAT-SDSS Galaxy Clusters Survey. 
Correlating X-ray and optical properties. 

P. Popesso, H. Bohringer, W. Voges 
Max-Planck-Institut fur extraterrestrische Physik, 85748 Garching, Germany 

Abstract. For a detailed comparison of the appearance of cluster of galaxies in 
X-rays and in the optical, we have compiled a comprehensive database of X-ray 
and optical properties of a sample of clusters based on the largest available X-ray 
and optical surveys: the ROSAT All Sky Survey (RASS) and the Sloan Digital Sky 
Survey (SDSS). The RASS-SDSS sample comprises 114 X-ray selected clusters. 
For each system we have uniformly determined the X-ray and optical properties 
. For a subsample of 53 clusters we have also compiled the temperature and the 
iron abundance from the literature. The total optical luminosity can be determined 
with a typical uncertainty of 20% with a result independent of the choice of local 
or global background subtraction. We searched for parameters which provide the 
best correlation between the X-ray luminosity and the optical properties and found 
that the z band luminosity determined within a cluster aperture of 0.5 Mpc h^" 
provides the best correlation with a scatter of about 60%. The resulting correlation 
of Lx and L op in the z and i bands shows a logarithmic slope of 0.38, a value not 
consistent with the assumption of a constant M/L. Consistency is found, however, 
for an increasing M/L with luminosity. This is confirmed by our analysis of the 
cluster mass to light ratio for a subsample of 53 clusters. We found a significant 
dependence of M/L on the cluster mass with a logarithmic slope ranging from 
0.27 in the i and r bands to 0.22 in the z band. 
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1. Introduction 

Cluster of galaxies are the largest well 
defined building blocks of our Universe. 
They form via gravitational collapse of cos- 
mic matter over a region of several mcga- 
parsecs. Cosmic baryons, which represent 
approximately 10-15% of the mass con- 
tent of the Universe, follow dynamically 
the dominant dark matter during the col- 
lapse. As a result of adiabatic compression 
and of shocks generated by supersonic mo- 



tions, a thin hot gas permeates the cluster 
gravitational potential. For a typical clus- 
ter mass of 10 14 M© the intracluster gas 
reaches a temperature of the order of 10 7 
keV and, thus, radiates optically thin ther- 
mal bremsstrahlung and line radiation in 
the X-ray band. In addition to the hot, 
diffuse component, baryons are also con- 
centrated in the individual galaxies within 
the cluster. These are best studied through 
photometric and spectroscopic optical sur- 
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veys, which provide essential information 
about luminosity, morphology, stellar pop- 
ulation and age. 

Since solid observational evidences indi- 
cate a strong interaction between the two 
baryonic components, one of our goals is 
the comparison of the X-ray and the opti- 
cal appearance of the clusters. We want in 
particular find optical parameters that pro- 
vide the closest correlation to the X-ray pa- 
rameters, such that we can predict within 
narrow uncertainty limits the X-ray lumi- 
nosity from these optical parameters and 
viceversa. On the basis of these consider- 
ations, we have created a large database 
of clusters of galaxies based on the largest 
available X-ray and optical surveys: the 
ROSAT All Sky Survey (RASS), and the 
Sloan Digital Sky Survey (SDSS). By care- 
fully combining the data of the two surveys 
we have created the RASS-SDSS galaxy 
cluster catalog. The X-ray-selected galaxy 
clusters cover a wide range of masses, from 
groups of 10 125 Mq to massive clusters of 
10 15 M G in a redshift range from 0.002 to 
0.45. The RASS-SDSS sample comprises all 
the X-ray selected objects already observed 
in the sky region covered by the Sloan 
Digital Sky Survey. For all derived quan- 
tities, we have used Ho = 70 kms _1 Mpc _ , 
O ro = 0.3 and n x = 0.7. 

2. The data. 

In order to create a homogeneous catalog 
of X-ray cluster properties, we have calcu- 
lated all X-ray parameters using only RASS 
data for all clusters in the sample. The 
X-ray luminosity has been calculated with 
the growth curve analysis (GCA) method 
used for REFLEX and NORAS 2, based 
on the RASS3 database (Bohringer et al., 
2001). The optical photometric data were 
taken from the SDSS (York et al. 2000 
and Stoughton et al. 2002). The SDSS con- 
sists of an imaging survey of ir steradians 
of the northern sky in the five passbands 
u, g, r ,i, z, in the entire optical range. 
Since the galaxies do not have sharp edges 
or a unique surface brightness profile, it is 



nontrivial to define a flux for each object. 
The SDSS photometric pipeline calculates 
a number of different magnitudes for each 
object: model magnitude, Petrosian mag- 
nitudes and PSF magnitudes. In the data 
analysis of this paper we used the Petrosian 
magnitudes for galaxies brighter then 20 
mag and the psf magnitudes for objects 
fainter than 20 mag. 

3. Optical Luminosity from SDSS 
data 

The total optical luminosity of a cluster has 
to be calculated after the subtraction of the 
foreground and background galaxy contam- 
ination. Since we have used only photomet- 
ric data from the SDSS galaxy catalog, we 
have no direct information on the cluster 
galaxy memberships. There are two differ- 
ent approaches to overcome this problem: 
a color cut or a statistical background sub- 
traction. We chose the latter approach since 
the former method may introduce a bias 
against bluer cluster members observed to 
have varying fractions due to the Butcher- 
Oemler effect. 

We have considered two different ap- 
proaches to the statistical subtraction of 
the galaxy background. First we have cal- 
culated a local background within an an- 
nulus centered on the cluster X-ray center. 
The annulus has been then divided in 20 
sectors and those featuring a larger than 3<r 
deviation from the median galaxy density 
are discarded from the further calculation 
(fig.lU left panel). 

As a second method we have derived a 
global background correction. The galaxy 
number counts N^ g (m)dm was derived 
from the mean of the magnitude number 
counts determined in five different SDSS 
sky regions, each with an area of 30 deg 2 
(figffl right panel). The source of uncer- 
tainty in this second case is systematic and 
originates the presence of large-scale clus- 
tering within the galaxy sample, while the 
Poissonian error of the galaxy counts is 
small due to the large area involved. After 
the background subtraction we found that 
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Fig. 1. Left panel: The plot shows a cluster region with the local background. The dots 
represent the galaxies in the sample. The biggest dots correspond to the brightest galax- 
ies in apparent r magnitude. The local background number counts have been calculated 
inside the annulus with inner radius equal to r2oo + 0.2 deg and a width of 0.5 degree. 
The regions with voids due to lack of data or close clusters have to be discarded in the 
background estimation. Right panel: Global background number counts as a function of 
magnitude in the r band. The error bars include the contribution of large scale struc- 
ture. The line shows the counts-magnitude relation expected for a homogeneous galaxy 
distribution in a universe with euclidian geometry: N(r) = A r 10°' 6 ( r ~ 16 ). The value of 
A r = 5.99 ± 0.52(0.5mag)~ 1 deg~ 2 is the results of the fit in Yasuda et al. 2001 for 
12 < r < 17. 



the signal to noise in the u band was too 
low to be useful, and performed our anal- 
ysis on the 4 remaining Sloan photometric 
bands g,r,i,z. 

In order to calculate the total cluster lu- 
minosity, we have calculated first the abso- 
lute magnitude, corrected for Galactic ex- 
tinction and K-correction. We have esti- 
mated L by using the (background cor- 
rected) magnitude number counts of the 
cluster galaxies with the following prescrip- 
tion: 

L = ^2 Ni(m)li(m) + / <f){m)dm (1) 

The sum on right side is performed over all 
the N magnitude bins with galaxy number 
iVj(m) and mean luminosity k(m). The in- 
tegral is an incompleteness correction due 



to the completeness limit of the galaxy 
sample at mu m = 21 mag in the five Sloan 
photometric bands. (j)(m) is the individual 
Schechter luminosity function fitted to the 
galaxy sample of each cluster. The incom- 
pleteness correction is of the order of 5-10% 
in the whole cluster sample. This means 
that the galaxies below the magnitude limit 
do not give a significant contribution to the 
total optical luminosity. Therefore the most 
important source of error is due to the con- 
tribution of the background galaxy num- 
ber counts. The uncertainty in each bin 
of magnitude is given by the Poissonian 
error of the bin counts (y/Ni ot (m), with 
i = 1,...,N) and the background subtrac- 
tion in each magnitude bin (05 (771)). Since 
the galaxy counts in the bins are indepen- 
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dent, the error in the luminosity is given 
by: 

N 

AL = (Y / (Nt ot (™)+<J l bg (mf)^. (2) 



4. Correlating X-ray and optical 
properties 

For a cluster in which mass traces optical 
light (M/L op t is constant), the gas is in 
hydrostatic equilibrium (T oc AT 2 / 3 ), and 
Lx oc T 3 , we expect the X-ray bolometric 
luminosity to be related to the optical lumi- 
nosity as L op oc L x 5 and to the intracluster 
medium temperature as L op oc T x ^. 

In this section we show that tight corre- 
lations exist between the total optical clus- 
ter luminosity and the X-ray cluster prop- 
erties as the X-ray luminosity and the intr- 
acluster medium temperature. 

To search for the best correlation be- 
tween optical and X-ray properties and to 
optimally predict for example the X-ray lu- 
minosity from the optical appearance, we 
are interested in an optical characteristic, 
which shows a minimum scatter in the X- 
ray/optical correlation. Therefore, we per- 
formed a correlation using 4 of the 5 SDSS 
optical band, g, r, i and z, to find out 
which filter should be used in the predic- 
tion. We used a fixed aperture to calculate 
the optical luminosities for all the clusters, 
to make no a priori assumption about the 
cluster size. Moreover, to check whether the 
scatter in the correlation depends on the 
cluster aperture, we did the same analy- 
sis with different apertures. To quantify the 
L op —Lx and the L op —T x relations, a linear 
regression in log-log space was performed. 
In any photometric band the scatter has a 
clear dependence on the cluster aperture by 
showing a region of minimum on the very 
center of the cluster, around 0.5 Mpc hy . 

Fig. shows the L op — Lx and L op — T x 
relation respectively, at the radius of mini- 
mum scatter, 0.5 Mpc h^ 1 , in the z band. 
In fact the i and the z bands have a slightly 
smaller scatter than in the other optical 



bands at any radius. The best fit in the z 
band for the L op — Lx and the L op — T x 
relations at the radius of minimum scatter 
are respectively: 

T IT i n ll-79±0.02 t / \0.45±0.03/qn 

Lop/^Q — 1U ^X\Rosat) (o) 

L op /L Q = 10 11 - 75±0 - 02 ix( Boi ) - 38±0 - 02 (4) 

L op /L Q = I0ll-42±0.06 T ^.12±0.08 (g) 

5. Discussion and conclusions 

The value of the exponent in the power law 
for the L opt — Lx(boi) relation, is around 
0.38 in the region of minimum. The val- 
ues are not consistent within the errors 
with the value of 0.5 predicted under the 
assumption of hydrostatic equilibrium and 
constant mass to light ratio. The same con- 
clusion can be derived for the L opt — Tx 
relation. A simple reason of the disagree- 
ment could be due to the assumption of a 
constant mass to light ratio. In fact, Girardi 
et al. (2002) analysed in detail the mass to 
light ratio in the B band of a sample of 
294 clusters and groups , finding M/L oc 
£0.33±o.03 rp ne same results was found by 

Lin et al. (2003) in the K band. Thus, if we 
consider this dependence of M/L from the 
optical luminosity with the assumptions of 
hydrostatic equilibrium, the new expected 
relation between the optical luminosity and 
the X-ray luminosity and temperature are 
L op oc L X A and L op oc T x 25 , respectively, 
which are in good agreement with our re- 
sults. The behavior of the M/L as a func- 
tion of the cluster mass and luminosity is 
confirmed also by our data as shown in 
fig. OH We used the subsample of cluster 
with known temperature to calculate the 
mass within rsoo with the M — T relation 
of Finoguenov et al. (2001). M/L shows a 
clear dependence on the cluster mass with 
a slope from 0.27 in the i and r bands to 
0.22 in the z band. 

The scatter in the L op — L x relation for 
the aperture with the best correlation (0.5 
Mpc hyQ 1 ), in the L op variable is 0.20, and 
the scatter in the Lx variable is 0.22 in the 
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Fig. 2. Left panel: Correlation between optical luminosities and X-ray luminosities. The 
fit is performed with a linear regression in the log(L opt ) — log(Lx) space for each of the 
4 optical bands. The solid and the dashed lines arc the results of the orthogonal and 
bisector method respectively. The error bars are at the 68% confidence level in both the 
variables. Right panel: the same as in the left panel but for the correlation between Lx 
the the ICM temperature T\ 



correlations obtained in i and the z bands. 
Therefore, by calculating the total cluster 
luminosity in the central part of the sys- 
tem, one can use the i or z band to predict 
the X-ray luminosity from the optical data 
with a mean error of 60%. In the same way 
the optical luminosity can be derived from 
Lx with the same uncertainty. Analogous 
results are obtained for the L op — T x rela- 
tion. Since the observational uncertainties 
in the optical and in the X-ray luminosity 
are about 20%, the scatters of 60% in both 
relations should be intrinsic. 

To better understand the relation be- 
tween the X-ray and the optical properties 
and their physical implications, the optical 
luminosity has to be calculated within the 
physical size of the systems. This will be 
discussed in a forthcoming publication. 
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Fig. 3. Left panel: cluster mass within rsoo as a function of the optical luminosity 
calculated within the same radius. Right panel: cluster mass to light ratio as a function 
of the cluster mass. Obvious consequence of the M — L op relation is that M/L shows 
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0.22 in the z band. This is in agreement with the prediction based on our new Lx — L op 
relation. 
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